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© The phase of received waves is detected at every symbol period T by means of a phase detector (12) on 
^ the basis of the phase of a local signal and the detected phase is inputted to a serial circuit of delay circuits 
^ (15i-15 N ) each of which lags a signal by a delay time T. Then signals with phases * n (n=0, 1,..., N) delayed by 

1 to N symbols are outputted and supplied to a metric calculating section (16). By adding the sum of a partial 
<v> sequence (A<t>i\ i = n + 1-q, n + 2-q,..., n) out of the candidates (A4> n ; n = 0, 1, .... N) of an N-symbol phase- 
00 difference sequence to the phase * n - q detected at a point of time (n-q)T (q = 1, 2, .... N), to determine the 
^ estimated value of the phase * n of the received waves. The v-th power of the absolute value of the difference 
|s* u n (q) between the estimated value and the phase of the received waves is used as the metric of the q-symbol 
CO phase-difference detection. E|u n (q)| v = \ n , the sum of the metric from q = 1 to q = n is used as a branch metric at 
© a point of time nT for the candidates of the phase-difference sequence. A pass metric A = EX n are found by 

calculating the sum of the branch metrics from a point of time 1T to a point of time NT for all the candidates of 
[Jj the N-symbol phase-difference sequence, and the sequence candidates which give the minimum value is used 

as a decoding sequence. 
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Technical Field 

The present invention relates to a differential detecting method for differentially detecting a digital signal 
transmitted as a phase difference sequence in a symbol interval and for obtaining a decoded sequence, in 
5 particular, to a differential detecting method and a differential detector thereof corresponding to a maximum 
likelihood sequence estimation technique. 

Related Art 

w Phase modulated waves are conventionally demodulated by coherent detection and differential detec- 
tion. In the coherent detection, the receive side reproduces a carrier wave as a reference signal, measures 
the phase of a received wave corresponding to the reference signal, and estimates a transmitted code. In 
this case, since the absolute phase is unknown, the sender side generally uses differential phase shift- 
keying modulation (DPSK) that modulates information corresponding to the variation of the phase of the 

75 carrier wave. Since the reproduced reference signal is not affected by noise and the like, low error rate can 
be accomplished. 

On the other hand, as the differential detection, differential phase detection and quadrature differential 
detection have been widely used. In the differential detection, the reference wave is formed of a received 
wave with a delay of one symbol interval. Thus, since no carrier wave reproducing circuit is required, the 

20 detecting circuit can be simply constructed and the detecting operation can be performed at high speed. 
Consequently, the differential detection is suitable for receiving a burst signal in time division multiple 
access (TDMA) communication. However, since the signal with a delay of one symbol interval is used as 
the reference signal, the reference signal tends to be adversely affected by thermal noise and the like. 
Thus, the error rate of the differential detection degrades in comparison with that of the coherent detection. 

25 As a result, depending on whether the detecting circuit is complicated, the burst signal is received, and so 
forth, one of the coherent detection and the differential detection is selected. 

For example, for four-phase DPSK, at a bit error rate 0.1%, the difference in bit energy-to-noise rate 
(Eb/No) between differential detection and coherent detection is 1.8 dB. To reduce the difference, a 
maximum likelihood quadrature differential detection for estimating a transmitted data sequence has been 

30 proposed as in Reference 1. The quadrature differential detector is composed of delay devices and 
multipliers. (Reference 1: D. Divsalar and M. K. Simon, "Multiple-symbol differential detection of MPSK," 
IEEE Trans. Commun., vol. 38, pp. 300-308, March 1990.) In addition, a technique for recursive estimation 
using the Viterbi algorithm has been proposed as in Reference 2. (Reference 2: D. Makrakis and K. Feher, 
"Optimal noncoherent detection of PSK signals," Electronics Letters, vol. 26, pp. 398-400, March 1990.) 

35 Assuming that an N-symbol phase difference sequence A</> n (where n = 1, 2, .... N) is being 
transmitted, the received signal is quadrature-differentially detected and a maximum likelihood sequence 
estimation is applied. M-phase DPSK signals that are received in an interval (n-1)T ^ t < nT can be given in 
the complex representation as follows. 

40 z(t) = (2Es/T) ,/2 exp j [<f> n + Q] + w(t) (01 ) 

where 4> n = {2nWM; m = 0, 1 M-1} is a modulated phase; Es is energy per symbol; T is one symbol 

interval; e is a phase difference between a received wave and a locally oscillated wave of the receiver; w(t) 
is noise of the receiver; and A<t> n = <f> n - 4>n-i is the n-th phase difference. After z(t) is filtered, it is sampled 

45 in the symbol interval. The obtained signal sample sequence is denoted by {Z n ; n = 0, 1 N}. In the 

technique of Reference 1, a sequence that maximizes a metric given by the following equation is selected. 

A = |z N + z N _i exp j A<f> N + z N - 2 expj(A</> N + A^ N -,) + ... + zo exp j (A<f> N + A^-i + ... + A<Jn)| 2 

(02) 

50 

Equation (02) can be modified as the following equation. 
N n 

A_= Z Re[z n {2 z n . q exp j (A<t>„ ♦ AO„-i ♦ --- + 

n=l q=l (03) 
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where Re[.] is the real part of a complex number; and ( )* Is a conjugate complex In equation (03) 
assuming that the upper limit of the summation with respect to q is L (< N) and the following equafon ,s 
defined as a branch metric 

L 

A n = Re[Z„{2 z n .„ exp j (A<t> B + + + ^ n -,-i))*l 

(04) 



70 



l L_1 states 



75 



20 



phase difference sequences are successively estimated corresponding to Viterbi algorithm with M 
as the technique of Reference 2. «t on 

In other words, in the maximum likelihood sequence estimation according to Reference 1, metrics of an 
N-sequences of <N-M) received signal samples in the symbol interval T are calculated and a sequence ^wrth 
the maximum metric is output. Thus, the number of times of the metric calculation becomes M . In other 
words, the number of times of the calculation at each time is M M /N. 

On the other hand, in the maximum likelihood sequence estimation correspond.ng to he Viterb. 
decoding described in Reference 2. there are M L "' phase difference sequences. A maximum kelihood I path 
at each time is selected. Since each state has arriving paths from M states out of M states at the jus 
preceding time, the number of times of the metric calculation becomes M x M = M . Thus the number 
of times of the metric calculation does not depend on the length of the transmitted symbol sequence 
Consequently, the number of times of the metric calculation of Reference 2 is much reduced in comparison 
with that of Reference 1. However, as the modulation level M of the modulation increases, the amount of 
25 calculating process exponentially increases. ^ «„^„*,oi 

An object of the present invention is to provide a maximum likelihood decodmg and different 
detecting method and a differential detector thereof that provide an error rate equal to or lower than a 
conventional maximum likelihood sequence estimation method with lesser calculation complexity jo Jhat 
much lower error rate with the same calculation complexity as the conventional maximum likelihood 
sequence estimation method without a tradeoff of maximum transmission data rate. 



30 



35 



40 



50 



55 



Disclosure of the Invention 

A first aspect of the present invention is a differential phase detecting method of an M-phase DPSK 
modulated signal, comprising the steps of detecting a phase * n of a received wave in ^P redete " n,n ^ 
transmitted symbol interval T corresponding to a local signal at a time nT. where n is any integer, tracing 

back a candidate of N-symbol phase difference sequences {A* n ; n - 1. 2 N} for q symbols so as to 

form a partial sequence {A* i; i = n, n-1, .... n + 1-q} and adding the sum of the partial sequence {A^i = 
n n-1 .. n-M-q} to a detected phase *„-, of q symbols before so as to obtain an estimated value *„• of 
the phase defining the v-th power value of the absolute value of a difference u n (q) between the 
estimated value *, and the phase *„ as a metric of a q-symbol differential phase detect,™ where v is a 
real number that is 1 or greater, adding the metric from q = 1 to n so as to obtain the o.tow,ng branch 
metric X„ = |u„d)T + IM2)| V + - + adding the branch metric from n-1 to N so as o obtain a 

path metric A = X, + K 2 + ... + x N for the candidate of phase difference sequences {A<(,„, i - 
I and defining an N-symbol phase difference sequence with the minimum path metric as a decoded 
sequence and outputting the decoded sequence. 

A second aspect of the present invention is a differential phase detecting method of an M-phase DPSK 
modulated signal, the differential phase detecting method using a path memory and a path metric memory, 
the path memory being adapted for storing both states (where Q is a predetermined in teger ^ hat ,s 2 
or greater) defined by Q modulated phase differences at each time and surviving M° paths each 
representing most likelihood path to each of the M°- states, the path metric memory be .ngj ^aped for 
storing a path metric that represents likelihood of a sequence to each state, the d.fferent.al P^^etecUng 
method comprising the steps of detecting a phase * n of a received signal in a predetermined transmitted 
symbol interval T corresponding to a local signal at a time nT, where n is any integer, tracing back a state 
S , of the M°-' states at a time (n-1)T along a surviving path stored in the path memory for an (L-Q) time, 

obtaining a sequence {A4> n -,; i = 1, 2 L-1} along the surviving path to the last state of a state , and 

adding a last symbol of a phase difference A^ at the time nT to the sequence so as to form a candidate 
sequence {A*„ ,- i = 0. 1, 2 L-1}. where L is a predetermined integer and L i Q. adding a detected 
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phase ¥ n _ q at a time (n-q)T to the sum of the phase differences of a partial sequence {A<f> n -i; i = 0, 1, .... 
q-1} of the candidate sequence so as to obtain an estimated value of the phase ^ n and calculating the 
difference between the estimated value and the phase tf n so as t0 obtain a phase error u n (q). adding the v- 
th power value of the absolute value of the phase error u n (q) from q = 1 to L so as to obtain the following 
5 branch metric 

X (S„-,-S n ) = + |Un(2)| v + + K(L)| V 

that represents likelihood of M branches from M states out of the M 0 " 1 states at the time (n-1)T to one 

io of the M 0 " 1 state S n at the time nT. The M branch metrics X(S n -i — Sn) to each state at the time nT are 
added to path metrics A(S n -i) of the state S„-i at the time (n-1)T that are read from the metric memory so 
as to obtain path metrics A(S n | S n -i) of the M candidate sequences that pass through the M states of S n -i, 
and compare the path metrics A(S n | S„-,) so as to obtain a state S n -i* with the minimum value. This state 
S n -i' represents a state at the time (n-1)T of a most likelihood path (i.e. a surviving path) to the state S n at 

75 the time nT. In this manner, path metrics of surviving paths for all the M Q-1 states at the time nT are 
obtained and compared to obtain a state S n ' with the minimum value. The path memory is traced back to 
the state S n - D for a predetermined interval DT starting from the state S„\ outputting, as a decoded symbol, 
a phase difference A<j> n - 0 that is one of Q-1 phase differences which construct the state S n - D - 

A third aspect of the present invention is a differential phase detecting method of an M-phase DPSK 

20 modulated signal, comprising the steps of detecting a phase * n of a received wave in a predetermined 
transmitted symbol interval T corresponding to a local signal, obtaining the difference * n - * n - q between the 
detected phase * n and a phase * n _ q of up to L symbols before, where q = 1, 2, .... L, obtaining the sum 
5 n _t(q) = X A^ n _i (where E is the sum of the phase differences from i = tto q-1) of the phase difference 
determined for each of up to (q-1) symbols before and obtaining the v-th power value of the absolute value 

25 or the difference of the difference u n (q) between the detected phase difference * n - * n - q and the sum of 
the added value 5 n - q (q) and the candidate phase difference A<*> n . as a metric of the q-symbol differential 
phase detection, adding metrics of the L phase differences so as to obtain a branch metric Xn = |u n (l)| v + 
+ |u n (L)| v for the candidate phase difference A</> n ', and outputting the phase difference alternative with 
the minimum branch metric as a determined phase difference A<fr n . 
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Brief Description of Drawings 



Fig. 1 is a schematic diagram showing state transition of a differential detecting method for explaining a 
differential phase detecting method according to a first embodiment of the present invention; 
35 Fig. 2 is a block diagram showing the differential phase detector according to the first embodiment; 

Fig. 3 is a schematic diagram showing state transition of a differential phase detecting method according 
to a second embodiment of the present invention; 

Fig. 4 is a graph showing a bit error rate in a computer simulation according to the first embodiment; 
Fig. 5 is a graph showing a bit error rate in a computer simulation according to the second embodiment; 
40 Fig. 6 is a schematic diagram showing state transition for explaining a differential phase detecting 
method according to a third embodiment of the present invention; 

Fig. 7 is a graph showing a bit error rate in a computer simulation according to the third embodiment; 
Fig. 8 is a block diagram showing a quadrature differential detector according to a fourth embodiment of 
the present invention; 

45 Fig. 9 is a table for comparing computation complexities corresponding to the third and fourth 
embodiments of the present invention and to References 1 and 2 in the case that sequence length N is 
L; 

Fig. 10 is a block diagram showing a differential phase detector according to a fifth embodiment of the 
present invention; 

so Fig. 1 1 is a graph showing bit error rate in a computer simulation according to the fifth embodiment; 

Fig. 12 is a block diagram showing a quadrature differential detector according to a sixth embodiment of 
the present invention; 

Fig. 13 is a graph showing bit error rate in a computer simulation according to the sixth embodiment; 
and 

55 Fig. 14 is a graph showing bit error rate in a computer simulation according to a seventh embodiment of 
the present invention. 
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Best Modes for Carrying out the Invention 

First, embodiments of a differential phase detecting method according to the present invention will be 



described. 



w 



(1) First Embodiment 

It is assumed that the carrier phase at time n-1 modulated by M-phase DPSK modulation method is 
denoted by 4> n -v The phase difference representing information to be transmitted at time n »s denoted by 
A<f> n = 2nWM (where m = 0, 1, .... M-1). The carrier phase * n at time n is given by 4>„-i + jW>r, An N- 

symbol phase difference sequence {A* n ;n = 0, 1 N} is transmitted. When an M-phase DPSK signal .s 

received in an interval (n-1)T £ t < nT, it can be expressed as follows. 



z(t) = (2Es/T) 1/2 expj (<*>„ + *) + w(t) 



(1) 



where = 2mWM (where m = 0. 1 M-1) is a modulated carrier phase; Es is energy per symbol; T is 

one symbol length; e is the phase difference between the received wave and the output of the loca 
oscillator of the receiver; w(t) is noise of the receiver; and A* n = * n -* n -i is the n-th phase difference of 
information transmitted. 

After the received signal z(t) is filtered and the outband noise is removed, a phase difference between 
the received signal z(t) and the local signal of the local oscillator is detected. The phase difference at time n 
is denoted by * n . The phase difference * n is given by the following equation. 

*n = <*>n + v« + 0 mod 2* (2) 

where „„ is phase noise caused by thermal noise; mod 2* is modulo operation, in which (x + 2tt) mod 2* = 

X (W Now as^shown in Fig. 1, it is assumed that a phase sequence of (N + 1) phases *o to * N are obtained 
at times 0 to N, respectively. The transmitted phases of symbols transmitted from the transmitter side are 
30 expressed by the following equations. 



20 



25 



35 



40 



4> N = 0 H - 1 + A 4> N 

<f> N- I = 4> N-t + A 4> N- I 

0 N - Q +, = <f> N-*+ A <t> W-qfl 



(3) 



(3) 



<f> x = <f> o+ A 4> i 

45 From equations (3), the following equations are satisfied. 

A A*n-i + — + A*h-«+i mod 2tt 

50 0n-0n-«= " 1: * A * »- ■ mod 2tt (4) 

i=0 

On the receiver side, the variation A* N (q) of the phases of the received signals in an interval from t = (N- 
55 q)T to NT can be given by the following equations. 
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A¥ N (q)=1 r N-1 r H^ mod 2k 

q-i 

s = £ A 0N-i + A ?n(q) mod 2k (5) 

i-0 

where A^q) = i> N - w-q is the phase difference noise. At this point, the fixed phase difference d between 
io the received signal and the local oscillator of the receiver is removed. Since i? N and 7j N - q can be 
approximated by independent Gaussian noise of average value = 0 and variance = 2a 2 , Aif N (q) becomes 
gaussian noise with average value = 0 and variance = 2a 2 . When a phase error vector is defined as u = 
{U N (1 ), u N (2) Un(N)} t (where T is a transposed matrix), the phase noise Arj N is expressed as follows: 

75 q-1 

/ZN (q)=A Y„(q)- L A#N-i mod 2k (6) 

i=0 

20 A joint probability density function P of u under the condition that the phase difference sequence A<j> = - 
(A<fri , A<fc A£ N ) T is transmitted can be given by the following equation. 



25 



30 



1 T R ~ 1 A* 

PC A* i A* )= exp { } 

(2^) N/8 (detR) 1 ^ 2 



(7) 



where R is N x N covariance matrix of Aij = {At> n (1), Ai? n (2), At 7n (N)} t ; and det R and R" 1 are 
determinant and inverse of a matrix R, respectively. When M phase differences take place with the same 
probability (normally this condition is satisfied), the maximum likelihood estimation of the transmitted phase 
35 difference sequences is performed by finding the phase difference sequence A<f>' = (A<m\ Ate', A^n') 1 " 
that maximizes the probability of equation (7). This is equivalent to find the phase difference sequence A<f>' 
that minimizes u T R~ 1 u as in equation (8). 



40 A0' = M I N ^ t R-*ac (8) 

over A <f> 

To find A<f>\ it is necessary to obtain R. The phase difference noise is Gaussian noise with the following 
45 properties. 

(AtjnO)) = 0 (9) 
<A7/ N (i)A7, N (j)> = 2a 2 for i = j 
= a 2 for i * j 

50 

where <x(i)> is the statistical expression of x. Thus, R can be given by the following equation. 
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R = (n j)= a 



15 



2 1 
1 2 



I 1 



1 2 1 
- 1 2J 



(10) 



When the standard matrix theory is applied, FT 1 can be given by the following equation. 
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r n -i 

l N 



R-«= <r -«/(N+l) 



•1 N -1 
- -1 NJ 



(11) 



When R"' is substituted into equation (8), the maximum likelihood estimation is performed by finding a 
phase difference sequence with the minimum value of the following equation from M N phase difference 
35 sequence candidates as shown in Fig. 1 . 



40 



N n 

A= £ £ I i.(q) I * 

n=l q=l 



(12) 



A computer simulation shows that even if |n„(q)| v (where v is a real number that is 1 or greater) instead of 
|u n (q)p is used, the similar result can be obtained. Thus, the present embodiment is arranged to find a 
45 phase difference sequence which minimizes the following value. 



50 



N n 

A= £ £ 



I a< »( q ) I 



(13) 



Fig. 2 shows a construction of a differential phase detector according to the first embodiment of the 
present invention. 

55 A received signal is supplied from an input terminal 11. The phase difference between the received 
signal and the local signal of a local oscillator 13 is detected by a phase detector 12. In the case of a 
transmitting device for use with present mobile communication, the received phase modulated signal is an 
intermediate frequency signal having a center frequency at the second intermediate frequency, and is the 



BNSOOCID: <EP 0671837A1J_> 



EP 0 671 837 A1 



output signal from a limiter amplifier or an AGC amplifier. The detected output is sampled by a sampling 
circuit 14 in a predetermined interval (symbol interval T). The sampled signal is input to delay circuits 15i 
to 15 N each of which has a delay of the symbol interval T. Thus, phases * n (where n = N, N-1, .... 0) with 
delays of 0 to N symbols are input to a metric calculating portion 16. 

s In the first embodiment, whenever a predetermined number N of phases of the received signal z(t) are 
received in the symbol interval T, the maximum likelihood estimation is performed for a sequence of N 
phase differences A£ n (where n = 1 to N). In other words, the metric calculating portion 16 calculates a 
path metric of phase difference sequence candidate corresponding to equation (12) and stores the path 
metric in a metric memory 16A. The metric calculating portions 16 repeats this process for all the 

w candidates (namely, 4 N in the case of four-phase DPSK), obtains a phase difference sequence with the 
minimum path metric, and outputs the result as a decoded output to a terminal 17. In reality, the metric 
calculating portion 16 adds the sum of partial sequence {A<fo; i = n + 1-q, n+2-q, n} of the N-symbol 
phase difference sequence candidate {A<J>„; n = 1, 2, .... N} to a detected phase * n - q at a time (n - q)T 
(where q = 1, 2, n) and obtains an estimated value V of a received signal phase * n corresponding to 

/5 the following equation. 

= *n-q + A*n + Mn-i + ••• + A<f> n - q+1 mod 2tt (14) 

A difference (estimated error) u n (q) between the estimated value * n ' and the received signal phase * n 
20 is obtained by the following equations from q = 1 to n corresponding to equation (6). 

25 ^,(2) = ¥n-(Yn-*+A* 0 + Atf B -i) mod 2k 

M n (q) = y B - % + & <f> n+ & 4> & & + mod 2tt 

/<,(n)=¥.-(*i+A^+A^-.+- + A*i) mod 2tt 
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(15) 

Next, it is assumed that the v-th power value of the absolute value of each estimated error u n (q) is a 
metric of q-symbol differential phase detection. The summation of the metrics from q = 1 to n at a time nT 
is given by the following equation. 



q=l 



As the most important point of the present invention, the phase difference at any time n is estimated 
not only from the phase *„_i of the just preceding signal, but also from phases * n -2. •■ of further preceding 
signals, and thus the bit error rate can be improved using these estimated errors u n (q). As in equation (12), 
so v is theoretically 2. However, simulation results show that v can be any positive rea! number in the range 
from 1 to 10. By the summation of the branch metrics from time 1T to time NT, the path metric A of the 
phase difference sequence candidates A<Jm, A<j> 2 , .... A<f> N is obtained corresponding to the following 
equation. 



55 
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A- L A. 07) 
n=l 

Path metrics are calculated for all (M N ) N-symbol phase difference sequence candidates and then 
stored in a metric memory 16A. A phase difference sequence candidate with the minimum path metric is 
defined as a decoded sequence and output from a terminal 17. 

(2) Second Embodiment 

In the first embodiment, the error rate are improved in proportion to the value of N. However, in the M- 
phase DPSK system, since there are M phase points, path metrics for all the M N sequence candidates 
shown in Fig. 1 should be calculated. Thus, the computation complexity is exponentially proportional to the 
value of N. In a second embodiment of the present invention, to satisfy the improvement of the error rate 
and the reduction of the computation complexity, the upper limit of the addition with respect to q in 
equation (13) is set to L (where L is any integer; 2 * L < N). In addition, the following equation is defined as 
a branch metric. 

L 

An = £ I mM \ ' (18) 
q=l 

As shown in a state transition of Fig. 3, phase difference sequences are successively estimated using 
Viterbi algorithm with M Q ~ 1 states. In this case, states S„ at a time nT are defined by a phase difference 

sequence {A4 n ,A4> n ~, A* n - Q+2 } from time n-Q + 2 to time nT. Since each phase difference has M 

values, the number of states at each time is M Q_1 . The relation between Q and L is 2 < Q £ L. 

The construction of a phase difference detector according to the second embodiment will be described 
with reference to Fig. 2. In this embodiment, the number of delay circuits 15 connected in series is L rather 
than N. Next, the process algorithm performed in the metric calculating portion 16 shown in Fig. 2 will be 
described. 

Step SI: When metrics X„(S„-i - S n ) of M branches that from M states S n -i = (A<fr n - (Q -i), 
A4> n - 2 , A* n -,) at a time (n-1)T to a state S n = (A<f> n - (Q - 2) , A* n -i, A* n ) at a time nT are calculated, 
surviving path to the state Sn-! is traced back for a (L-Q) time and a sequence of (L - 1) symbols (A*„-i; i 

= 1,2 L-1) is obtained. A phase difference A<f> n at a time nT that is a last symbol is added to the 

symbol' sequence to obtain an L-symbol sequence and then the metric X n <S n _i - S n ) is calculated 
corresponding to equation (18). The number of states at each time is M Q ~\ However, the number of states 
S,,-! at the time (n-1)T that can be shifted to each state S n at the immediately succeeding time nT is M out 
of M 0 " 1 (where M is the number of phases of the DPSK system). 

Step 2: A path metric of M branches to the state S n is calculated corresponding to the following 
equation. 

A n (S n | S n -i) = A n -,(S n -,) + X n (S„-i-S n ) (™) 

The path with the minimum path metric is a surviving path with the maximum likelihood that arrives at the 
state S n ,and the state from which the transition has been made is stored in the path memory 16B (see Fig. 
so 2). The value of the minimum path metric is defined as a path metric x n (S n ) of each state and stored in the 
metric memory 16A. 

Step S3: A path that gives the minimum path metric A n (S n ) among M Q states S n is selected as a 
surviving path and the state transition stored in the path memory 16B is traced back along the selected path 
for an interval DT. The depth of the trace back may be about 4. Phase difference A<j> n - D of state 
55 (A4> n - <D+Q -2). • A* n _ (D+1) , A* n - D ) traced back are output as decoded results A</> n - D -. 

Fig. 4 shows the results of a computer simulation of the error rate characteristics of the differential 
phase detection according to the first embodiment. In the simulation, v = 2. In the graph, the horizontal axis 
represents bit energy-to-noise ratio (Eb/No). In Fig. 4, the cases of N = 2, 3, and 4 are plotted by O , A, 
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and respectively. In addition, for comparison, error rates of conventional symbol-by-symbol differential 
phase detection (N = 1) and coherent detecting with differential decoding are plotted by x and +, 
respectively. In Fig. 4, solid lines 21 and 22 represent ideal curves. To accomplish the error rate of 0.1%, 
the difference in Eb/No between symbol-by-symbol differential phase detection and the coherent detecting 

5 with differential decoding is 1 .8 dB. When N = 3, the difference in Eb/No can be reduced to 0.9 or less. 
When N = 4, the difference in Eb/No becomes around 0.6 dB. In Fig. 4, the theoretical characteristics 
(Reference 1) of error rate in the case that the maximum likelihood sequence estimation is performed for 
the conventional quadrature differentially detected output are denoted by dotted lines. From Fig. 4, it is 
clear that in the first embodiment, almost the same improvement as Reference 1 can be accomplished. 

w Fig. 5 shows the results of a computer simulation of the error rate characteristics of the differential 
phase detection method according to the second embodiment. When L = Q = 4, the difference in Eb/No 
between the symbol-by-symbol differential phase detection and the coherent detecting with differential 
decoding becomes approximately 0.2 dB. 

In the second embodiment, the relation between Q and L is 2 £ Q £ L. Next, a third embodiment where 

75 Q = 2 and Q < L will be described. 

(3) Third Embodiment 

In the second embodiment, the Viterbi algorithm with M Q_1 states (where 2 £ Q £ L) was used so as to 
20 reduce the computational complexity in comparison with that of the first embodiment. In the third 
embodiment, a differential phase detection method with lesser computational complexity will be described. 
As with the second embodiment, in the third embodiment, the maximum likelihood sequence estimation is 
performed based on the Viterbi algorithm. As shown in a state transition diagram (M = 4) of Fig. 6, in the 
third embodiment, the number of states of the Viterbi algorithm is M that is the same as the number of 
25 phases of modulation. In other words, in the third embodiment, Q = 2. Thus, a state Sn at a time nT is Atf> n . 
The phase difference represents the state. Consequently, there are only M survival paths at each time. In 
the second embodiment, the number of survival paths at each time is M°~\ Next, the construction of the 
third embodiment will be described with reference to Fig. 2. 

First, a phase * n of a received signal z(t) with respect to a locally oscillated signal is detected by the 
30 phase detector 12 in the transmitted symbol interval T. As described above, a detected phase at a time nT 
is * 0 = <*>n + 0 + Vn where ij n is phase noise caused by thermal noise. (L + 1) received phase samples 

{* n - q ; q = 0, 1 L} are used for the predetermined value L. As is clear from equation (5), * n and * n _ q 

have the following relation. 

35 

Y. 1" Arf-i + C*.-*.-,) mod 2k (20) 

i-0 

40 

Phase errors u n (1). Un(2) a n (L) between the received phase * n and the estimated phase * n * are 

calculated corresponding to the following equation. 



*-(q )=¥.-¥.-,- 1 A0»-i mod2;r (21) 

i=0 

As with equation (18), the sum of the v-th power values of the absolute values of the phase errors are 
50 calculated from q = 1 to L as a branch metric \„. Thus, the branch metric can be given by the following 
equation. 

L 

A (A 0 A 4> »)= 2 U.(q)l' (22) 

q=l 
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where v is any real number that is 1 or greater. With the branch metric expressed by equation (22), 
symbols are decoded using Viterbi algorithm with M states (that will be described later). 

Step S1 : To select a most likely path that arrives at a state A*„ at a time nT from M phase difference 

states at a time (n-1)T, a phase difference sequence {A*„-,; I - 1. 2 L-1) stored in the path memory 

16B is read by tracing back surviving paths for a past time (n-L + 1) starting from one state A4>„-, of wi 
states at a time (n-1)T. The state A<f>„ at the time nT that is a last state is added to the phase difference 
sequence. Thus, the phase difference sequence candidate {A<K-,; i = 0. 1 L-1} is formed. 

Step S2: A detected phase *„-„ at a time (n-q)T is added to the sum of the phase differences of the 

partial sequence {A*,,-* i = 0. 1 q-1} of the phase difference sequence candidate so as to obtain an 

estimated value *„' of a phase *„. A phase error u.(q) that is the difference between the estimated value 
* • and the detected phase *„ is obtained from q = 1 to L corresponding to equation (15). 

Step S3- The v-th powers of the absolute values of the L phase errors u„(q) are calculated. The 
resultant values are added from q = 1 to L corresponding to equation (22). Thus, a branch metric MA*„-, 
- A</>„) that represents likelihood of transition from the state A*„_, at the time (n-1)T to the state A</>„ at tne 
time nT is obtained. Next, as in the following equation, the branch metric is added to the path metnc A- 
(A*„-,) of the state A<f>„-, at the time (n-1)T stored in the metric memory 16A so as to obtain a path metric 
A(A4>„ | A<f>„-,) of a candidate sequence that passes through the state A<f>„-1 



A(A* n | A* n -,) = A(A4>„-,) + \(A<f>„~i-A<j> n ) 



-1- 

(23) 



Step S4- The above-described calculations are performed for each of the M states A<f> n -, at the time (n- 
1)T so as to obtain path metrics for the M candidate sequences leading to the state A*„. By comparing the 
path metrics, a state A*„-,' with the minimum value is obtained. The state A4>„-i' that is decided as a state 
at the time (n-1)T on the most likely path to the state A* n at the time nT is stored in the path memory 16B. 
25 The path metric A(A*„ | A*„-,') is decided as a path metric A(A* n ) of the state A*„ at the time nT and is 
stored in the metric memory 16A. ...... A . „♦ t u D 

Step S5- The processes and calculations at the steps SI to S4 are repeated for all M states A4>„ at the 
time nT so as to obtain M path metrics. By comparing the M path metrics, a state A*„' with the minimum 
value is obtained. The path memory is traced back for a predetermined interval DT starting from the state 
30 A<*>„' and the obtained state A<f>„- D is output as a decoded symbol. 

Fig. 6 shows an example of a state transition diagram in the case of M = 4 according to the third 
embodiment. In this example, the maximum likelihood path to the state A*„ = 0 at the tome nT is selected 
Referring to Fig. 6, paths denoted by dotted lines extend from four states A*,-, = 0, W2, and at 
the time (n-1)T to A<f» n = 0. Survival paths denoted by solid lines extend to each state at the time (n-i)T. 
35 For example, in the calculation for the path metric including a transition branch from one state A<J>„-, - »/Z 
at the time (n-1)T to one state A* „ = 0 at the time nT. the state is traced back for (L-1) states along , a 
survival path (stored in the path memory 16B) that extends to the state A<f>„-, = W2 at the timeji,-^ 
Thereafter, a branch metric MA*„-, - A*„) is calculated corresponding to equations (15) and (22). The 
branch metric is added to the path metric A(A<j>„-i) of the state A* n -, = »/2 at the time (n-1)T stored in the 
40 metric memory 16A so as to obtain a path metric A(A<j>„ | A* n -,) of the candidate sequence. Such a 
process is repeated for the M states A« n -, at the time (n-1)T so as to obtain respective path metrics. By 
comparing the M path metrics, the maximum likelihood path to the state A<o n at the time nT ,s selected 

Fig 7 is a graph showing the results of a computer simulation of error rate performances of me four- 
phase DPSK method according to the third embodiment. In this simulation, v = 2. In the graph ^the 
horizontal axis represents the ratio of signal energy and noise power density per bit (Eb/No) and the vertica 
axis represents the error rate. For comparison, in Fig. 7, simulation results of error rates of the conventional 
symbol-by-symbol differential phase detection (L = 1) and the coherent detecting with differential decoding 
are plotted. In the graph, solid lines represent theoretical curves. The difference in Eb/No between he 
symbol-by-symbol differential phase detection and the coherent detecting with differential decoding at the 
error rate of 0.1% is 1.8 dB. However, when L = 2, the difference in Eb/No is the half or less of 1.8 dB. In 
Fig 7 theoretical performance curves of the error rate of the maximum likelihood sequence estimation of 
Reference 1 are plotted with dotted lines. When L = 2 and 4, performances equivalent to L - 3 and 5 in 
Reference 1 are accomplished, respectively. As the value of L increases, the performance close to 

Reference 1 is accomplished. . »i 

Thus in the case of performing the most likelihood sequence estimation through the V.terb. algonthm, 
when the number of states of the Viterbi decoder is equal to the number of phase of the modulation, the 
computational complexity can be remarkably reduced in comparison with that of Reference 1. 
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(4) Fourth Embodiment 

In the third embodiment, the phase difference between the received signal z(t) and the locally oscillated 
signal of the local oscillator 13 is detected by the phase detector 12 so as to perform the Viterbi decoding 
5 with the M states. However, as with the third embodiment, the Viterbi decoding with M states can be 
performed for a sample sequence of a complex detection output of a quasi-coherent detection for the 
received signal z(t). This method will be described as a fourth embodiment of the present invention. Fig. 8 
is a block diagram showing the construction of a quadrature differential detector according to the fourth 
embodiment. 

w A received signal z(t) is supplied from an input terminal 11 to a quasi-coherent detecting circuit 12. The 
quasi-coherent detecting circuit 12 performs quasi-coherent detection for the input signal with two locally 
oscillated signals supplied from a local oscillator 13. The frequency of the input signal is nearly the same as 
that of the locally oscillated signals. The locally oscillated signals are 90 deg. apart in phase. The complex 
output of the quasi-coherent detecting circuit is supplied to a sampling circuit 14. The sampling circuit 14 

75 samples the complex signal in a predetermined interval (symbol interval T) and outputs a complex sample 
Z n of the received signal. The complex sample Z n is input to delay circuits 15i to 15 L that are connected in 
series and each of which causes delay of a symbol period T. The delay circuits 15i to 15 L output samples 

with delays of 1 to L symbols {Z n _ q ; q = 1, 2 L}. The delayed samples and non-delayed samples Z n 

are input to a metric calculating portion 16. The metric calculating portion 15 has a metric memory 16A and 

20 a path memory 16B. The metric calculating portion 16 performs calculations corresponding to a decoding 
algorithm similar to that in the third embodiment. The decoding algorithm will be described in the following. 
Decoded output data is obtained from a terminal 17. 

Step S1: To select a most likely one of the paths arriving at a state A<f» n at a time nT from M phase 
difference states at a time (n-1)T, a surviving path stored in the path memory 16B is traced back to a past 

25 time (n-L+1) starting from a state A<fr n -i of M states at the time (n-1)T so as to read a phase difference 

sequence {A*,,-*; i = 1, 2 L-1}. The state A<j> n at the time nT as the last state is added to the phase 

difference sequence so as to form a phase difference sequence candidate {A4> n -i; i = 0, 1 L-1}. 

Step S2: The phase of the received signal sample Z n - q is rotated for the sum of the partial sequence 
{A<f> n _j; i = 0, 1 q-1} of the sequence candidate. This process is repeated from q = 1 to L The L 

30 obtained values are added so as to obtain an estimated value z n ' of the received signal sample z n . 

Step S3: As expressed in equation (04), a real value of the inner product of the received signal sample 
z n and the estimated value z n ' is defined as a branch metric M&tn-} — A</> n ) that represents the likelihood 
of the transition from the state A<t> n --\ at the time (n-1)T to the state A£ n at the time nT. The branch metric 
is added to a path metric ti(A<t> n -i) of the state A^ n _, at the time (n-1)T as in equation (23) so as to obtain 

35 a path metric A(A</> n | A^n-,) of a candidate sequence that passes through the state A^n-,. 

Step S4: The above-described calculations are repeated for the M states A^-t at the time (n-1)T so as 
to obtain path metrics for the M candidate sequences. By comparing the path metrics, a state A<t> n -i' with 
the maximum value is obtained. The state A</> n -r is decided as a state at the time (n-1)T of the most likely 
path to the state A<f> n at the time nT and stored in the path memory 16B. The path metric A(A<j> n | A<j> n -0 is 

40 defined as a path metric A(A<*> n ) of the state A£ ft at the time nT and stored in the metric memory 16A. 

Step S5: The processes and calculations at the steps S1 to S4 are repeated for all M states A* n at the 
time nT so as to obtain M path metrics. By comparing the M path metrics, the state A<f>„' with the maximum 
value is obtained. The path memory is traced back for a predetermined interval DT starting from the state 
A^r/ and the obtained state A<*> n - 0 is output as a decoded symbol. 

45 In the quadrature differential detecting based on the Viterbi decoding described in Reference 2, there 
are M L_1 survival paths at each time. However, in the methods according to the third and fourth 
embodiments of the present invention, since there are only M survival paths, the error rate thereof are 
slightly inferior to that of Reference 2. However, in these embodiments, the calculating amount of the 
branch metrics at each time is only M 2 . Thus, the computational complexity in the third and fourth 

50 embodiments is much smaller than that of Reference 2. Fig. 9 shows the number of calculations of the 
branch metrics at each time in the case of M = 4. According to the results of the simulation of the method 
proposed in Reference 2, when the Viterbi algorithm with M L ~ 1 states is used, the error rate performance 
equivalent to 2L in Reference 1 can be accomplished. The number of calculations of the branch metrics for 
obtaining the error rate performance equivalent to the case of L = 6 in Reference 1 is compared among the 

55 methods of References 1 and 2 and the present invention. The number of branch metric calculations in 
Reference 1 is 683; the number of calculations in Reference 2 is 4096; and the number of calculations in 
the present invention is 16. Thus, it is clear that the number of calculations in the present invention can be 
remarkably reduced in comparison with that of the conventional methods. 
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(5) Fifth Embodiment 

In the third and fourth embodiments, the Viterbi algorithm is used. Thus, one most likely path leading to 
each of a predetermined number of states at each time from the just preceding time is selected. Symbol 
corresponding to states at a time traced back for a predetermined interval along the most hkely path 
se ected amo'ng them are output as a decoded result. Thus, the decoded result has a delay agains a 
received phase sample for the predetermined number of symbols. As a maximum take hood decoding 
method that is free of such a delay, a decision feedback decoding algonthm that does not use the Jiterbi 
algorithm is known. In the decision feedback decoding algorithm, there is only one suiting pa* at each 
time. Based on the single path of a sequence that has been determined the "^JT^" 
determined. The decision result is immediately output. Next, the dec.s.on feedback decoding method will be 
described as a fifth embodiment of the present invention. .... .,• t 

Fig 10 is a block diagram showing the construction of a differential phase detecting c.rcurt according to 
the fifth embodiment of the present invention. A phase modulated signal z(t) is ; received frorr .an input 
terminal 11. The phase modulated signal z(t) is supplied to a phase detector 12. The phase «J ofto 
phase modulated signal z(t) is detected corresponding to the phase of a locally oscillated signal of a local 
oscillator 13. An output of the phase detector 12 is supplied to a sampling circuit 14 The samphng wcurt 
14 samples the input signal in a symbol interval T and outputs a phase sample *„ o th ^ e '^ signal to 
a differential phase detecting portion 15. In the differential phase detecting portion 15 the detected phase 

*„ of each symbol is supplied to L delay circuits 15,. 1* 15 L each of which has a delay of one symbol 

interval T. The differences *„ - ♦„-,. *„ - *« - *n- L between the delay outputs *„-,, *„- a . . 

of the delay circuits 15,, 15, 15 L and the input detected phase *„ are obtained by sub racting circuits 

15S ,. 15S 2 15S L . The phase differences accord with those expressed by equation (5). The phase 

differences are supplied to a metric calculating portion 16. 

The metric calculating portion 16 decides on the phase difference A*„ corresponding to a calculating 
process as will be described later and outputs it to an output terminal 17. In addition, the phase difference 
M n is supplied to a cumulating portion 18. In the cumulating portion 18. the detected phase difference A*„ 

is input to (L-1) delay circuits 18,, 18* 18,-,. each of which is connected in series and has a delay of T. 

The delay circuits 18,, 182. .... 18 L -, supply delay outputs A4„-„ A*„- 2 t0 add '" 9 A C '™ 

18A, 18A 2 I8A3 , 18A L - 2 . respectively. Outputs of the adding circuits 18A,, 18A2 18A l - 2 are 

successively supplied to next-stage adding circuits 18A 2 , I8A3 18A L _ 2 . Outputs of the .add. ng circurts 

18A, to 18A,- 2 are supplied to the metric calculating portion 16. In other words, $„-,(q) - EA*„-, (where l 
is the sum from i = 1 to q and q = 1 to L-1) is supplied to the metric calculating portion 16 The metric 
calculating portion 16 selects one candidate from M phase differences A*»' (in the case of four-phase 
DPSK system, 0, W2, », and 3*/2) and adds the selected one to each added value 5„-,(q) supplied from 
the cumulating portion 18. The v-th power value of the absolute value or the difference u„(q) (corresponding 
to equation (6)) of the selected candidate phase difference alternative and the detected phase difference 
(*n-* q ) is obtained as follows. 

|U„(1)| V - |*n-*n-1 -A*n'| V (24) 
|M-„(2)| V = |*„-*„-2-(A*„'jrA<J. n -,2mod27r| v 
|U.o(3)| V = |* n -+ n -3-(A*n'_+A<».n-, + A^^-2) mod 2ir| v 

|U„(L)| V = l+^n-L-tA^' + A*,,-,*-*- +A<<. n - L+ ,)mod27r| v 

where v is a real number that is 1 or greater. The sum K = qu„(q)r of |u„d)r to |u n (L)r is defined as the 
branch metric of the candidate phase difference A*„\ For all M phase difference candidates ^n. their 
branch metrics are calculated and the candidate phase difference A« n ' with the minimum branch metric is 
output as the detected phase difference A<f.„. 

Fiq 11 is a graph showing results of a computer simulation of error rate performance of four-phase 
so DPSK system according to the fifth embodiment. In this simulation, v = 1. In the graph, the horizontal axis 
represents bit energy-to-noise ratio (Eb/No). In Fig. 11, the performance in the case of L - 1 is the same 
as that of the conventional symbol-by-symbol phase delay detection. In Fig. 11, for comparison, the enw 
rate performance of the coherent detecting and differential decoding is shown. The difference ,n Eb/No 
between the symbol-by-symbol differential phase detection and the coherent detecting with differential 
55 decoding at the error rate of 0.1% is 1.8 dB. When L = 3. the difference in Eb/No can be reduced to 
almost the half of 1.8 dB. When L = 10, the difference in Eb/No becomes 0.2 dB. 
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(6) Sixth and Seventh Embodiments 

In the third and fourth embodiments, the number of states of the Viterbi decoding algorithm is reduced 
to M (in other words, the number of calculations of the branch metrics is. reduced). In addition, states are 

5 traced back for L-2 along surviving paths leading to each state A</> n -i at the time (n-1)T so as to obtain 
phase difference sequence from A^„_ 2 to A<j> n _ L +i. Branch metrics are calculated using the state A^ n at the 
time nT and past (L-1) states and a maximum likelihood sequence estimation is applied based on the 
Viterbi algorithm. Thus, the error rate is improved. In sixth and seventh embodiments according to the 
present invention, branch metrics can be calculated corresponding to a plurality of past phase states without 

w need to perform the trace-back process at each time. Thus, in the sixth and seventh embodiments, the 
number of calculations can be more reduced than that of the third and fourth embodiments. 

Equation (02), which expresses a path metric in the maximum likelihood sequence estimation for the 
above-described quadrature differential detection is modified as follows. 

75 N n 

A = £ Re[z»{E z„-<exp j (A 4> + A 4> n- q+ i ))*exp- j A 4> » J 
n=l q=l 

(25) 

20 

The portion in brackets { } of equation (25) is defined as a reference signal h,,^ for calculating a branch 
metric. The reference signal h n _ t is given by the following equation. 

25 n 

hn-i= L z„.,exp j (A <f> +A <f> „-, + ■) (26) 

Q-l 

30 Thus, the branch metric X of equation (04) can be given by the following equation. 

X = RelZnhViexp-jA^] (27) 

According to equation (26), the reference signal h^ is the sum of samples from time (n-1)T to 0. When a 
35 forgetting factor 0 (0 ^ 0 £ 1) is used in equation (26) so as to make the contribution of past samples to 
hn-T to be in proportion to time, the following equation is obtained. 

40 h.-,= Z j^-'z— ,expj (A*.-i+-+A *.-t + i) (28) 

q=i 

The phase difference partial sequence in the bracket ( ) of equation (28) is traced back from time n-1. Thus, 
45 the phase difference in the case of q - 1 is 0. Consequently, in the region of q > n + 1 (that is, before 
time 0), if z n - q = 0, equation (28) can be given by the following recursive expression. 

oo 

h B -i = z„-i + E $ «- | z.-,exp j +A <f> „-< + i) 

q=2 



50 



55 



= z„-i+ £h n - 2 expj A <f> n-« (29) 



In equation (27), the branch metric for the state transition from the state A<fr n -i to the state A<fr n is given by 
the following equation. 
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x(A*„-i-A*„) = Re[z n hVi(A*n-i)exp-jA4>„] (30) 

where h n -i(A* n -i) is the value of h„-, for a state A4> n _i at a time n-1. Using a reference signal h n - 2 (A4>„-. 2 ) 
for a state A* n -2 at a time n-2 of a surviving path to a state A4>„-i at a time n-1, a reference signal h„-r 
(A^ n -i) for a state A<fr n -i can be calculated in the recursive manner. Thus, unlike with the third and fourth 
embodiments, in the sixth and seventh embodiments, it is not necessary to trace back a survival path to the 
state A^n-i so as to calculate the reference signal h,,-, corresponding to equation (26). In addition, in the 
sixth and seventh embodiments, the delay circuit used in the third and fourth embodiments is omitted. Next, 
a decoding algorithm corresponding to the sixth and seventh embodiments will be described. 

Fig. 12 is a block diagram showing the construction of a quadrature differential detector according to 
the sixth embodiment of the present invention. As with the embodiment shown in Fig. 8, a received signal z- 
(t) is quadrature-detected by a quasi-coherent detecting circuit 12. The quasi-coherent detecting circuit 12 
outputs complex detection output to a sampling circuit 14. The sampling circuit 14 samples the detector 
output every symbol interval T. The sampling circuit 14 outputs a received signal complex sample z n . In the 
sixth embodiment, the delay circuits 15, to 15 N in Fig. 8 are not used. The received signal complex sample 
z n is supplied to a metric calculating portion 16 that comprises a metric memory 16A, a path memory 16B, 
a branch metric calculating portion 16C, a Viterbi algorithm portion 16D, a reference signal calculating 
portion 16E, and a reference signal memory 16F. The metric calculating portion 16 estimates a transmitted 
phase difference sequence based on the Viterbi algorithm with M states in the following steps. 

Step 1: To select a most likely path to a state A<fr n at a time nT from M phase difference states at a 
time (n-1)T, a value h(A4> n -,> of a reference signal h^ of a state A^„-, of M states at the time (n-1)T is 
read from the reference signal memory 16F. The branch metric calculating portion 16C calculates a branch 
metric \(A<f> n -i - A<f> n ) that represents likelihood of transition from the state A<fr n -t at the time (n-1)T to the 
state A4> n at the time nT corresponding to equation (30). 

Step S2: A path metric A(A4> n -i) of the state A<f> n -i at the time (n-1)T is read from the metric memory 
16A. The Viterbi algorithm portion 16D adds the branch metric X(A<j> n -i - A<f> n ) to the path metric A(A<f> n -i) 
of the state A*„-, at the time (n-1)T so as to obtain a path metric A(A<f>„ | A$ n -i) of a candidate sequence 
that passes through the state A<£ n -i. 

Step S3: The processes at the steps 1 to 3 are repeated for all M states A* n -i so as to obtain path 
30 metrics of the M candidate sequences. The Viterbi algorithm portion 16D compares path metrics of the M 
candidate sequences and obtains a state A<fr„-r with the maximum value. The Viterbi algorithm portion 16D 
decides the state A^-V as a state of a surviving path at the time (n-1)T that leads to the state A* n at the 
time nT and stores it in the path memory 16B. The Viterbi algorithm portion 16D defines the path metric A- 
(A<f> n | A4>„-D as a path metric A(A<> n ) of the state A</> n at the time nT and stores it in the metric memory 
35 16A. 

Step S4: The reference signal calculating portion 16E calculates a value h n (A<f> n ) of a reference signal n n 
of a state A<f>„ used in the calculation at the next time (n + 1)T corresponding to the following equation and 
stores it in the reference signal memory 16E corresponding to the path metric A(A<f> n ). 

h n (A<f> n ) = Zn + jShn-^A^n-^expjA^n ( 31 ) 
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Step S5: The processes at the steps 1 to S4 are repeated for all the M states at the time nT so as to 
obtain M surviving paths and path metrics. By comparing the path metrics, a state A<#> n ' with the maximum 
value is obtained. The path memory 16 is traced back for a predetermined interval D starting from the state 
45 A<t> n \ The obtained state is output as a decoded symbol A</> n -D- 

Fig. 13 is a graph showing results of a computer simulation of error rate performance of a four-phase 
DPSK system according to the sixth embodiment of the present invention. In the graph, the horizontal axis 
represents the ratio of signal energy and noise power density per bit (Eb/No). For comparison, results of 
computer simulation for error rate of the conventional symbol-by-symbol differential detection and coherent 
detecting with differential decoding are plotted in Fig. 13. The difference in Eb/No between the conventional 
symbol-by-symbol differential detection and the coherent detecting with differential decoding at error rate of 
0.1% is 1.8 dB. However, when /3 = 0.9, the difference in Eb/No can be reduced to 0.1 dB or less. Thus, 
according to the sixth embodiment, the computational complexity can be more reduced than that of the 
third and fourth embodiments. When the forgetting factor fi is nearly 1, the error rate performance close to 
the coherently detecting and differential decoding method can be almost accomplished without need to 
increase the calculating amount. 

In the sixth embodiment, the number of survival paths at each time is M. However, in the seventh 
embodiment, the number of survival paths at each time is limited to one. Thus, the construction of the 
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seventh embodiment is simplified. The construction of a differential detector according to the seventh 
embodiment is basically the same as the construction shown in Fig. 12. However, in the seventh 
embodiment, the metric memory 16A and the path memory 16B are not required. A decoding algorithm of a 
metric calculating portion 16 according to the seventh embodiment is described in the following steps. 

5 Step SI: To decide which path arriving at one of M states A$ n at a time nT from the phase difference 
state A^ n -i decided at a time (n-1)T provides the maximum likelihood, a value h(A$ n -i) of a reference 
signal h n _! corresponding to the decided state at the time (n-1)T is read from the reference signal 

memory 16F. The branch metric calculating portiorM6C calculates, using the reference value (A<f> n -i) and 
the received signal sample z n , a branch metric X(A<^„_ 1 -* A</> n ) that represents likelihood of transition from 

w the decided state A^„_i at the time (n-1)T to one of the M states A<» n at the time nT based on equation 
(30). 

Step S2: The process at the step S1 is repeated for ail the M stages A<f> n so as to obtain branch 
metrics for the M candidate states. The Viterbi algorithm portion 16D compares the branch metrics and 
obtains a state A<> n with the maximum value and outputs the state A<f> n as a decoded symbol A<j> n . In this 
75 decoding method, the path metrics of all the states are 0. Thus, the symbol decision is performed by 
comparing the branch metrics. 

Step S3: The reference signal calculating portion 16E calculates a value h n (A$ n ) of a reference signal h„ 
of the state A^ n to be used for the calculation at the next time (n-1)T based on the following equation and 
stores the obtained value in the reference signal memory 16E. 
20 _ _ 

h n (A4> n ) = Zn + jShn-^n-^eXDjA^ (32) 

Fig. 14 is a graph showing results of a computer simulation of error rate characteristics of a four-phase 
DPSK system according to the seventh embodiment. In the graph, the horizontal axis represents the ratio of 

25 signal energy and noise power density per bit (Eb/No). As with the sixth embodiment, when 0 - 0.9. the 
difference in Eb/No between the conventional symbol-by-symbol differential detection and the coherent 
detecting with differential decoding can be reduced to 0.1 dB or less. 

As described above, in the differential detection method according to the present invention, the error 
rate can be much more improved than that of the conventional symbol-by-symbol differential detection. 

30 Thus, the error rate performance close to the coherent detection with differential decoding can be 
accomplished. Alternatively, when the error rate performance attainable with conventional method is to be 
attained, the calculating amount can be much reduced in comparison with the conventional method. In the 
first, second, third, and fifth embodiments of the present invention, since estimation of maximum likelihood 
phase difference sequence is performed using detected phase sequences of received signal, a practical 

35 limiter amplifier can be used. In addition, since a fast acquisition property that is a feature of the differential 
detection is not lost, the present invention can be applied to burst reception in the TDMA system. 

Claims 

40 1. A differential phase detection method of an M-phase DPSK modulated signal, comprising the steps of: 

a) detecting a phase <£ n of a received signal in a predetermined transmitted symbol interval T 
corresponding to a local signal at a time nT, where n is any integer; 

b) adding the sum of a q-symbol partial sequence {A<fr; i = n, n-1, .... n + 1-q} of an N-symbol 
phase difference sequence candidate {A<j> n ; n = 1, 2, N} to a detected phase * n - q of q symbols 

45 before so as to obtain an estimated value of the phase 

c) defining the v-th power value of the absolute value of a difference u n (q) between the estimated 
value * n * and the phase as a metric of a q-symbol differential phase detection, where v is a real 
number that is 1 or greater; 

d) adding the metric from q = 1 to n so as to obtain the following branch metric 

50 

x n = |M1)| V + |Un(2)| v + ... + |u„(n)| v ; 

e) adding the branch metric from n = 1 to N so as to obtain a path metric A = Xi + X 2 + ... + x N 
for the candidate phase difference sequences {A<f> n ; n = 1, 2, .... N}; and 

55 f) defining an N-symbol phase difference sequence with the minimum path metric as a decoded 

sequence and outputting the decoded sequence. 
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A differential phase detection method of an M-phase DPSK modulated signal, the differential phase 
detection method using a path memory and a path metric memory, the path memory being adapted for 
storing as a surviving path, both M Q_1 states, where Q is a predetermined integer that .s 2 or greater, 
defined by Q modulated phase differences at each time and states each representing from which one 
of states at just preceding time a most likely path arrives at each of the M° states, the path metric 
memory being adapted for storing a path metric that represents likelihood of a sequence leading to 
each state, the differential phase detection method comprising the steps of; 

a) detecting a phase of a received signal in a predetermined transmitted symbol interval T 
corresponding to a local signal at a time nT. where n is any integer; 

b) tracing back from one state S n -, of the M°-' states at a time (n-1)T along a surviving path stored 

in the path memory for an (L-Q) time, obtaining a sequence {A* n -,; i - 1,2 L-1} along the 

surviving path and adding, as a last symbol, a phase difference A* n at the time nT to the sequence 
so as to form a candidate sequence {A<f> n -,; i = 0, 1 , 2 L-1}. where L is a predeterm.ned integer 

padding a detected phase at a time (n-q)T to the sum of the phase differences of a partial 

sequence {A4> 0 - i: i = 0. 1 q-1} of the candidate sequence so as to obtain an estimated value of 

the phase and calculating the difference between the estimated value and the phase so as to 
obtain a phase error u„(q); 

d) adding the v-th power value of the absolute value of the phase error u„(q) from q - to L so as 
to obtain the M branch metrics that represent likelihood of M branches from the M«" states S„-, at 
the time (n-1)T to a state S„ at the time nT 

X n = |U„(1)| V + |U.„(2)| V + ... + |li n (n)| v ; 

e) adding the branch metric X(S n -, - S„) to a path metric AfS,.-,) of the state S„-, at the timejn-1)- 
T that is read from the metric memory so as to obtain path metric A(S„ | S n - ? ) of the candidate 
sequences that pass through the state S„-„ and comparing all of the M path metrics A(S„ | S„-,) so 
as to obtain a state S„-,' with the minimum value; 

f) defining the state 8^i" as a state at the time (n-1)T of a survival path leading to the state S n at the 
time nT, storing the state S„-,' in the path memory, defining the path metric A(S n | S„-,') as a path 
metric A(S„) of the state S n at the time nT, and storing the path metric A(S„) in the path metric 

gTrepelting the calculations at the above step for all the M Q_1 states at the time nT so as to obtain 
M Q -' path metrics and comparing the M 0-1 path metrics so as to obtain a state S„ with the 

35 minimum value; and , .... 

h) tracing back the path memory for a predetermined interval DT starting from the state S„ , defining 
a phase difference A^-d that is one of Q-1 phase differences that constructs the state S n - D as a 
decoded symbol, and outputting the decoded symbol. 

3. The differential phase detection method as set forth in claim 2, wherein L = Q. 

4. The differential phase detection method as set forth in claim 3, wherein Q = 2. 

5. The differential phase detection method as set forth in claim 2, wherein L > Q. 

6 A quadrature differential detection method of an M-phase DPSK modulated wave, the quadrature 
differential detection method using a path memory and a path metric memory, the path memory being 
adapted for storing both M states that represent modulated phase differences at each time and phase 
difference states at just preceding time of maximum likelihood paths each leading to one of the M 
states at each time, the path metric memory being adapted for storing path metrics that represent 
likelihood of sequences leading to each state, the quadrature differential detection method comprising 

^^sampling a received signal at a time nT in a predetermined transmitted symbol interval T so as to 
obtain a received signal sample z„; 

b) tracing back the path memory starting from one state A4>„-1 of the M states at a time (n-l)T, 
obtaining a sequence {A*,,-,; i = 1. 2. .... L-1} along a surviving path to the state A* n _,. add.ng a 
state A* n . and forming a candidate sequence {A<f>„-i; i - 0, 1. 2 L-1}; 
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c) rotating the phase of a received signal sample z n _ q for the sum of the partial sequence {A$„_j; i 
= 0, 1, 2, .... L-1} of the candidate sequence from q = 1 to L and adding L obtained values so as to 
obtain an estimated value z n ' of the received signal sample z n ; 

d) calculating a real value of the inner product of the received signal sample z n and the estimated 
value and defining the reaJ value as a branch metric X(A<f> n -i — A$ n ) that represents likelihood of 
transition from the state &<t> n -y at the time (n-1)T to the state A<f> n at the time nT; 

e) adding the branch metric \{A$ n -<i — A<*> n ) to a path metric A(A<^ n _ 1 ) of the state A^ n -! at the time 
<n-1)T so as to obtain a path metric A(A<f> n | A«> n -i) of the candidate sequence that passes through 
the state A<t> n -u 

f) repeating the above calculations for all the M states A^n-! so as to obtain path metrics of the M 
candidate sequences and comparing the M candidate sequences so as to obtain a state A<f> n _i' with 
the minimum value; 

g) deciding the state A<f>„-i' as a state at the time (n-1)T of a surviving path to the state A<fr n at the 
time nT. storing the state A<^ n -i' in the path memory, defining the path metric A(A<*> 0 | A^n-/) as a 
path metric A(A*„) of the state A<f> n at the time nT, and storing the path metric A(A<f> n ) in the path 
metric memory; 

h) repeating the calculations at the steps b) to g) for all the M states at the time nT so as to obtain M 
path metrics and comparing the M path metrics so as to obtain a state A<f> n ' with the minimum value; 
and 

i) tracing bac* the path memory for a predetermined interval DT starting from the state A4» n f and 
outputting an obtained state A</> n _ D as a decoded symbol. 

A differential phase detection method of an M-phase DPSK modulated signal, comprising the steps of: 

a) detecting a phase * n of a received signal in a predetermined transmitted symbol interval T 
corresponding to a local signal; 

b) obtaining the difference - *„_ q between the detected phase * n and a phase * n - q of up to L 
symbols before, where q = 1.2 J_; 

c) obtaining the sum & n -i(q) = E A<j> n - U where E is the sum of the phase differences from i = 1 to 
q-1 decided for each of up to (q-1) symbols before, and obtaining the absolute value or the v-th 
power value of the difference u n (q) between the detected phase difference * n - * n - q and the sum of 
the added value 5 n _ q (q) and the phase difference candidate A£ n ' as a metric of the q-symbol 
differential phase detection; 

d) adding metrics of the L phase differences so as to obtain a branch metric X n = |u n (1)| v + - + 
|u n (U| v for the phase difference candidate A<fr n '; and 

e) outputtiing the phase difference candidate with the minimum branch metric as a decided phase 
difference A<> n . 

A quadrature differential detecting method of an M-phase DPSK modulated signal, the quadrature 
differential detection method using a path memory, a path metric memory, and a reference signal 
memory, the path memory being adapted for storing both M states that represent modulated phase 
differences at each time and phase difference states at just preceding time of most likely paths each 
leading to one of the M states, the path metric memory being adapted for storing path metrics that 
represent likelihood of sequences each leading to one of M states, the reference signal memory being 
adapted for storing a reference signal used for calculating a branch metric for each state, the 
quadrature differential detection method comprising the steps of: 

a) sampling a received signal in a transmitted symbol interval T so as to obtain a received signal 
sample z n at a time n; 

b) phase-rotating a value hn-^A^n-i) of a reference signal \\ n -^ in one of M states A<f> n _i at a time 
(n-l)T by A<f> n> calculating a branch metric X(A4> n -, - A<f> n ) that represents likelihood of transition 
from the state A<f> n _, at the time (n-l)T to the state A«> n at the time n corresponding to a real value 
of the inner product of the phase-rotated reference signal value and the received wave sample z n , 
adding the branch metric to a path metric A(A<> n - 1 ) of the state &<t> n -, at the time (n-l)T, and 
obtaining a path metric A<A<f> n | A<f> n -i)of a candidate sequence that passes through the state A<f> n -i; 

c) repeating the process at the step b) for all M states A<f> n -i so as to obtain path metrics of the M 
candidate sequences; comparing the path metrics so as to obtain a state A4> n -i* with the maximum 
value, defining the state A<f> n -,' as a state at the time (n-1)T of a most likely path leading to the state 
A<f> n at the time n, storing the state A4> n -t* in the path memory, and defining the path metric A(A<*> n | 
A<J>n-0 as a path metric A(A<f> n ) of the state A$ n at the time n, and storing it in the metric memory; 
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d) calculating a value h n (A<f> n ) of a reference signal h n of the state A<f> n , the value h n (A<f> n ) being to be 
used for calculations at the next time (n+ 1)T, corresponding to the following equation 

h n (A* n ) = Z n + /Shn-^n-iJexpj A<frn 

5 

where /S is a predetermined constant and 0 < /3 £ 1 , and storing the value h n (A<f> n ) in the reference 
signal memory; and 

e) repeating the processes at the steps b), c), and d) for all M states at the time nT so as to obtain 
M surviving paths and path metrics, comparing the path metrics, obtaining a state A4> n ' with the 

w maximum value, tracing back the path memory for a predetermined time D starting from the state 

A<j> n \ and outputting the obtained state A<f> n - D as a decoded symbol. 

9. A quadrature differential detection method of an M-phase differential modulated signal, the quadrature 
differential detection method using a metric memory and a reference signal memory, the metric 
75 memory being adapted for storing a path metric that represents likelihood of a sequence leading to 
each state, the reference signal memory being adapted for storing a reference signal used for 
calculating a branch metric of each state, the quadrature differential detection method comprising the 
steps of: 

a) sampling a received signal in a transmitted symbol interval T so as to obtain a received signal 
20 sample z n at a time n; 

b) rotating the phase of a reference signal h n -! associated with a determined phase difference state 
A^ n _, at a time (n-1)T so as to_ determine a most likely one of M states A<fr n at the time n and 
calculating a branch metric \(A^ n - A -> A</> n ) that represents likelihood of transition from the state 
A^n-i at the time (n-1)T to the state A<fr n at the time n corresponding to a real value of the inner 

25 product of the rotated reference signal and the received signal sample z n ; 

c) repeating the above calculations for all the M states A<f> n so as to obtain branch metrics of the M 
sequences, comparing the branch metrics so as to obtain a state A$ n with the maximum value, 
determining the state A<f> n as a decoded jymbol A<f> n of the received signal sample z n at the time 
nT, and outputting the decoded symbol A<J> n ; and 

30 d) obtaining the value of a reference signal h n associated with the state A<f> n used for calculations at 

the next time (n + 1 )T corresponding to the following equation 

h n = z n + 0h n _i exp j A^ n 

35 where 0 is a predetermined constant and 0 < £ 1 , and storing the value in the reference signal 

memory. 

10. A differential phase detector, comprising: local signal generating means for generating a local signal 
with a frequency nearly the same as that of a transmitted signal modulated with M-phase DPSK at 
40 symbol intervals T; 

phase detecting means for receiving the transmitted signal and detecting the phase of the received 
signal corresponding to the phase of the local signal; 

sampling means for sampling the phase of the received signal every symbol interval T and 
outputting a sampled phase * n ; 
45 delay means having N delay stages connected in series and each causing delay for the symbol 

interval T, where N is an integer that is 2 or greater, said delay means being adapted for receiving the 
sampled phase from said sampling means and outputting sampled phases * n _i.*n-2, *n-N of 1 to N 
symbols before; and 

metric calculating means for receiving the present sampled phase from said sampling means and 
50 N past sampled phases from the delay stages of said delay means, adding the sum of a partial 
sequence {A*,; i = n, n-1, and n + 1-q} of q symbols before in an N-symbol phase difference 
sequence candidate {A<> n ; n = 1, 2, .... N} to a detected phase * n - q of q symbols before so as to 
obtain an estimated value * n * of the phase * n , defining the v-th power value of the absolute value of a 
difference u n (q) between the estimated value ¥ n ' and the phase ¥ n as a metric of a q-symbol phase 
55 detection, where v is any real number that is 1 or greater, adding the metric from q = 1 to n so as to 
obtain the following branch metric 

Xn = |Un(1)| V + M2)| V + + |Un(n)|\ 
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adding the branch metric from n = 1 to N so as to obtain a path metric A = X1 + X2 + ... + XN of 
the phase difference sequence candidate {A<f>j; i = 1, 2, .... N}, and outputting an N-symbol phase 
difference sequence with the minimum path metric as a decoded sequence. 

11. A differentia! phase detector for successively estimating a sequence of phase differences at symbol 
intervals of a phase signal obtained by phase-detecting an M-phase DPSK modulated signal based on 
a Viterbi algorithm with M Q_1 states defined by past (Q-1) phase difference sequence so as to decode 
the phase difference sequence, where Q is a predetermined integer that is 2 or greater, the differential 
phase detector comprising: 

local signal generating means for generating a local signal with a frequency nearly the same as that 
of a transmitted signal modulated corresponding to M-phase DPSK modulation system in a symbol 
interval T; 

phase detecting means for receiving the transmitted signal and detecting the phase of the received 
signal relative to the phase of the local signal; 

sampling means for sampling the phase of the received wave every symbol interval T and 
outputting a sampled phase *„; 

delay means having L delay stages connected in series and each causing delay for the symbol 
interval T, where L is an integer and Q £ L, said delay means being adapted for receiving the sampled 
phase from said sampling means and outputting sampled phases *„-,, * n -2. *h-l of 1 to L symbols 
before; and 

metric calculating means having a path memory and a path metric memory, the path memory 
being adapted for successively storing the states of survival paths at the just preceding time each 
leading to one of M Q_1 states, the path metric memory being adapted for storing path metrics that 
represent likelihood of surviving sequences each leading to one of M Q_1 states, said metric calculating 

means being adapted for receiving (L-1) sampled phases * n -i. *n-2 *n-L at past times (n-1)T, (n- 

2)T (n-L + 1)T from the delay portions of said delay means and decoding the sampled phases based 

on the Viterbi algorithm with the M Q_1 states. 

wherein said metric calculating means is adapted for adding the sum of a partial sequence {A</m ; i 
= n, n-1, .... n + 1-q} of q symbols of an L-symbol phase difference sequence candidate {A<f> n ; n = 1, 
2, L} to a detected phase *v q of q symbols before so as to obtain an estimated value * n ' of the 
phase * n , defining the v-th power value of the absolute value of a difference u n (q) between the 
estimated value V and the phase * n as a metric of q-symbol differential phase detection, where v is a 
real number that is 1 or greater, calculating a metric X n (S n -i - S n ) of M branches that leave M Q_1 
states Sn-, = (A<fr n _ (Q -t,, .... A^„- 2 , A<f> n -i> at the time (n-1)T to each state S n = (A<*> n - ( Q- 2 ). ».. A<f> n -i, 
A<f> n ) at the time nT according to the following equation 

X n = K0)| v + M2)| v + - + |Un(L)| v . 
calculating path metrics of the M branches to the states S n according to the following equation 

A n (S n | S n -i) = An-tfSn-,) + X^Sn-, - S„) 

defining a branch with the minimum path metric as a survival path to the state S n , storing the survival 
path in the path memory, defining the value of the minimum path metric as a path metric A n (S n ) of the 
state S n , storing the path metric A n (S n ) in the metric memory, selecting a path leading to a state with 
the minimum path metric A n (S n ) of all the M Q_1 states S n at the time nT, tracing back states stored in 
the path memory along the selected path for a DT interval, and outputting a phase difference A</> n -D of 
states (A4> n _ (D+Q -2), A<> n - (D+ i), A<^ n - D ) as a decoded result A^ n - D '- 

12. A differential phase detector, comprising: local signal generating means for generating a local signal 
with the same frequency as that of a received signal M-phase DPSK modulated at symbol intervals T; 

phase detecting means for detecting the phase of the received signal relative to the local signal 
and outputting the detected phase * n at the symbol intervals T; 

phase difference detecting means for obtaining a phase difference (* n -*n- q ) between the detected 
phase and a phase * n -q of up to L symbols before, where L is an integer that is 2 or greater; 

cumulating means for cumulating decided phase differences 5 n -i(q) = EA<J>n-i. where Z is the sum 
of the phase differences from i = 1 to q-1 ; 
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metric calculating means for calculating the v-th power value or the absolute value of the difference 
between the detected phase difference * n -¥„- q and the sum of the added value 5 n -i(q) and a phase 
difference candidate A* n ' as a metric X q of a q-symbol differential phase detection, adding L metrics of 
Ihe differential phase detection so as to obtain a branch metric X n = |u. n 0)| v + - + IMJ-)" and 
outputting a phase difference candidate with the minimum branch metric as a decided phase difference 
Ao«. 

13. The differential phase detector as set forth in claim 12, wherein said phase difference detecting means 
includes delay means having L delay stages and L subtracters, the L delay stages being connected in 
scr.es and each causing delay for a symbol interval T, the L subtracters being adapted for successively 
receiving detected phases from said phase detecting means and outputting past sampled phases *„-i, 

* n _ L of up to 1 to L symbols before, the subtracters being adapted for subtracting the output 

pnascs of the delay stages from the detected phase * n received from said phase detecting means and 
outputting the phase difference * n -*n- q - 

14. The different.al phase detector as set forth in claim 12, wherein said cumulating means has delay 
moans and <L-1) adders, the delay means having L delay stages connected in series and each causing 
delay lor a symbol interval T, the adders being adapted for receiving successively decided phase 
d ftofonces from said metric calculating means, receiving past decided phase differences A<f> n _i, A<*> n -2, 

20 A* n t of 1 to I symbols before from the delay stages, cumulating the output phases of the delay 

starts, and generating a sum 5 n -i(q) = EA*n-i from i = 1 to q-1 and from q = 2 to L. 

15. A quadrature differential detector for receiving and detecting a transmitted signal which has been M- 
phase DPSK modulated at symbol intervals T, comprising: 

local signal generating means for generating two local signals with nearly the same frequency as 
the transmitted signal and with a phase difference of 90 deg each other; 

quasi-coherent detecting means for receiving the transmitted signal, quadrature^ detecting the 
received signal with the two local signals, and outputting samples z n of a complex signal at symbol 
intervals T; 

a path memory for storing a phase difference state at the just preceding time of a most likely patn 
leading to each of M states that represent modulated phase differences at each time; 

a metric memory for storing a path metric that represents likelihood of a sequence to each state; 
a reference signal memory for storing a reference signal used for calculating a branch metric for 
each state* 

35 branch metric calculating means for receiving a received signal sample z n at a time n from said 

quasi-coherent detecting means, phase-rotating a value h^A*^) of a reference signal h n _! read out 
from said reference signal memory by A<f> n , and calculating a branch metric \(A</> n _i - A<fr n ) that 
represents likelihood of transition from the state A+„-, at the time (n-1)T to the state A</> n at the time n 
corresponding to a real value of the inner product of the rotated reference signal value and the received 

40 signal sample z„; 

Viterbi algorithm means for adding the branch metric to a path metric A(A4>„-i) of the state A<f>„-1 
at the time (n-1)T read from said metric memory, repeatedly generating a path metric A(A<*> n | A<f> n -i) of 
a candidate sequence that passes through the state A<j> n -, for all M states A*„-i. obtaining path 
metrics of all the M candidate sequences, comparing the path metrics, obtaining a state A*„-i' with the 

45 maximum value, deciding the state as a state at the time (n-1)T of a most likely path leading to 

the state A* n at the time n, storing the state in said path memory, defining the path metric A(A<f> n | 
A</> n -Oas a path metric A(A$ n ) of the state A<f> n at the time n, and storing the path metric A(A<f> n ) in 
said metric memory; 

reference signal calculating means for calculating a value h n (A<*> n ) of the reference signal h n 
50 associated with the state A<» n according to the following equation 

h n (A<f>„) = z n + j8h n -i(A4> 0 -i)expj A*„ 

where 0 < 0 S 1, and storing the value h n (A<f> n ) in said reference signal memory, 
55 wherein said Viterbi algorithm means is adapted for obtaining survival paths and path metrics for all 

the M states at the time nT, comparing the path metrics, obtaining a state A</> n ' with the maximum 
value, tracing back said path memory for a predetermined interval D starting from the state A*„\ and 
outputting an obtained state A<f> n -o as a decoded symbol. 
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16. A quadrature differential detector for receiving and detecting a transmitted signal which has been M- 
phase DPSK modulated at symbol intervals T, comprising: 

local signal generating means for generating two local signals with nearly the same frequency as 
the transmitted signal and with a phase difference of 90 deg each other; 
5 quasi-coherent detecting means for receiving the transmitted signal, quadrature^ detecting the 

received signal with the two local signals, and outputting samples z n of a complex signal at symbol 
intervals T; 

a path memory for storing a phase difference state at the just preceding time of a most likely path 
leading to each of M states that represent modulated phase differences at each time; 
jo a metric memory for storing a path metric that represents likelihood of a sequence to each state; 

a reference signal memory for storing a reference signal used for calculating a branch metric for 
each state; 

branch metric calculating means for receiving a received signal sample z n at a time n from the 
quasi-coherent detecting means, phase-rotating the value of a reference signal h n -, associated with a 

75 decided phase difference state A<f> n _j at a time (n-1)T read out from said reference memory by A4> n , 
and calculating a branch metric \{A4> n - A A<j> n ) that represents likelihood of transition from the state 
A*,,-, at the time (n-1)T to the state A<f>„ at the time n corresponding to a real value of the inner 
product of the rotated reference signal and the received signal sample z n ; 

Viterbi algorithm means for adding the branch metric XfA^n-, — A<f> n ) to the path metric A at the 

20 time (n-1 )T read from said metric memory, repeatedly generating path metrics A(A<f> n ) to the M states 
A</> n , comparing the path metrics, obtaining a state A<f> n with the maximum value, determining the state 
A<f> n as a decoded symbol A$„ of the receivedjsignal sample z n at the time nT, outputting the decoded 
symbol A*„, determining the path metric A(A<M as a path metric A, and storing the path metric A in 
said metric memory; and 

25 reference signal calculating means for calculating a value of the reference signal h n associated with 

the state A</> n to be used for calculations at the next time (n + 1)T corresponding to the following 
equation 

h n = Z n + 0h n -i expj A^ n 

30 

where 0 is a constant and 0 < j8 * 1 , and storing the value of the reference signal hn in said reference 
signal memory. 
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